EULER MACLAURIN WITH REMAINDER FOR A SIMPLE INTEGRAL 

POLYTOPE 



YAEL KARSHON, SHLOMO STERNBERG, AND JONATHAN WEITSMAN 

Abstract. We give an Euler Maclaurin formula with remainder for the sum of the values of a 
smooth function on the integral points in a simple integral polytope. We prove this formula by 
elementary methods. 



1. Introduction 

The Euler Maclaurin formula computes the sum of the values of a function / over the integer 
points in an interval in terms of the integral of / over variations of that interval. A version of this 
classical formula is this: 

For any function f(x) on the real line and any integers a < b, we will consider the weighted sum 

(1) E 'f ■= ^^(") + /(« + 1) + ••• + /(&- 1) + Ifib). 

[a,b] 

If / is "nice enough" , for instance, a polynomial, then 

(2) E'/ = 

[a,b] 

where 

(3) US) 

Because fl^^j!^^ f{x)dx is a polynomial in hi and /12 if / is a polynomial in x, applying the infinite 
order differential operator L(^) then yields a finite sum, so the right hand side of is well 
defined when / is a polynomial. 

A polytope in M" is called integral, or a lattice polytope, if its vertices are in the lattice ; it is 
called simple if exactly n edges emanate from each vertex; it is called regular if, additionally, the 
edges emanating from each vertex lie along lines which are generated by a Z-basis of the lattice Z". 

Khovanskii and Pukhlikov |KPHlKP2] . following Khovanskii |KhHlKh2] . generalized the classical 
Euler Maclaurin formula to give a formula for the sums of the values of polynomial or exponen- 
tial functions on the lattice points in higher dimensional convex polytopes A which are integral 
and regular. This formula was generalized to simple integral polytopes by Cappell and Shaneson 
inHHEHSinHSllHl, and subsequently by Guihemin [H^l and by Brion-Vergne [BVl. All of these 
generalizations involve "corrections" to the Khovanskii-Pukhlikov formula when the simple poly- 
tope is not regular. When applied to the constant function / = 1, these Euler Maclaurin formulas 
compute the number of lattice points in A in terms of the volumes of "dilations" of A. A small 
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sample of the literature on the problem of counting lattice points in convex polytopes is given in 
|pH iMdl IVl iKKl iMol IPoml IryR" "RPR R^; see the survey [BP| and references therein. 

These formulas are closely related to the Riemann Roch formula from algebraic geometry via the 
correspondence between polytopes and toric varieties. Under this correspondence, regular polytopes 
correspond to smooth toric varieties, and the Khovanskii-Pukhlikov formula was motivated in this 
way, although it was proved combinatorially. Cappell and Shaneson derived their formula from 
their theory of characteristic classes of singular algebraic varieties and had the key idea of using 
the operator L as we do here. Guillemin obtained his formula from the equivariant Kawasaki- 
Riemann-Roch formula and methods coming from symplectic geometry and the theory of geometric 
quantization. Brion and Vergne employed a method that is closer to that in the original proof of 
Khovanskii and Pukhlikov, using Fourier analysis. 

To illustrate the relation to toric varieties, let us sketch the symplectic- geometric proof for the 
case of a regular polytope, following Guillemin |Gulj . This approach will not be used elsewhere in 
this paper. A regular integral polytope A C M"" determines a smooth Kahler toric variety (M, w), 
and geometric quantization gives rise to a virtual representation Q{M) of the torus T". The 
dimension dim(5(M) of this quantization is equal to the number of lattice points in A. (This result 
(see [01 Corollary 2.23]) is an expression of the "quantization commutes with reduction" principle 
in symplectic geometry {(irSj . According to this principle, dimQ{MY = dim (5(Mc) for each lattice 
point c G C Lie(r"')*, where Q{MY is the subspace of Q{M) on which T" acts through the 
character given by c, and where is the reduced space of M at c. Because M is a toric variety, Mc 
is a point if c G A and is empty otherwise.) On the other hand, by the Hirzebruch-Atiyah-Singer 
generalization of the classical Riemann-Roch formula, we have d\m.Q{M) = exp(ci(L))Td(rM), 
where ci(L) = [uj] is the Chern class of the pre-quantization line bundle and Td(rM) is the Todd 
class of the tangent bundle. Expressing M as a reduction of a linear torus action on (where d is 
the number of facets of A), the tangent bundle stably splits into line bundles Li, . . . , L^, and the 
above integral is obtained by applying the Khovanskii-Pukhlikov differential operator nTd(^) 
to the integral J^jexp{uj + ^^jCi(Lj)). The Duistermaat-Heckman theorem on the variation of 
reduced symplectic structures implies that this integral is equal to the volume of the polytope A(/i) 
that is obtained from A by shifting the ith. facet by a distance hi, for i = 1,. . . ,d. Hence, the 
number of lattice points in A is obtained by applying the Khovanskii-Pukhlikov operator to the 
volume of A(/i). 

The Euler Maclaurin formulas due to Khovanskii-Pukhlikov, Cappell-Shaneson, Guillemin, and 
Brion- Vergne are all exact formulas, valid for sums of exponential or polynomial functions. Cappell 
and Shaneson |CS3j have also investigated the problem of deriving an Euler Maclaurin formula 
with remainder. In a previous paper |KSWj . we stated and proved an Euler Maclaurin formula 
with remainder for the sum of the values of an arbitrary smooth function on the lattice points in 
a regular polytope, and adumbrated a generalization to the case of simple integral polytopes. The 
purpose of this paper is to state and prove an Euler Maclaurin formula with remainder for simple 
lattice polytopes (Theorem |21). The key ingredients in the proof of this theorem, as in |KSWj . 
are a variant of the Euler Maclaurin formula in one dimension, given in Proposition \27\ which by 
iteration gives a formula for orthants, along with a combinatorial result, given in Proposition I38| 
which shows how the sum of the values of a function over the lattice points in a polytope can be 
decomposed into sums over such orthants. 

Our Euler Maclaurin formula with remainder is stated in Theorem |2l for functions of compact 
support. In Section [3 we show how to extend it to symbols (in the sense of Hormander, see e.g. 
|Hoj ) . This is the content of Theorem |21 As a corollary, we deduce an exact Euler Maclaurin 
formula for polynomials. 

The early references to the Euler Maclaurin formula are Euler |Euj and Maclaurin |Maj . Appar- 
ently, Poisson |Poij was the first to give a remainder formula. See also |Harj . 
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2. Weighted sums in one dimension 

Exact Euler Maclaurin. Here is a brief proof of the exact Euler Maclaurin formula ^ ; cf . 
First, we prove this formula when f{x) is an exponential function: f{x) = e^^ with |A| < 27r. The 
formula then becomes 

o o /'b+h2 

[a^b] -J a til hi=h2=0 

An explicit computation, which uses the facts that the constant term in the formal power series 
L(S) is one and that L(— 5) = L(S'), shows that 



(5) / e^^dx 



dhi dh2 



h+h 



a—h 



f-b 

=/l2=0 



where \-?'^{S) is the truncation of the power series L(S') at the even integer 2k. The radius of 
convergence of the power series L(A) is 2ti because the zeros of tanh(A/2) that are nearest to the 
origin are at zizlni. Hence, for |A| < 2tt, the expression in (jSJ converges as A; ^ oo to 

(0) L(A) /' .-dx = f^! _ = i f + e- " 



tanh(A/2) V A A; 2V ; _ g- V2 ' 

which is equal to the left hand side of Q by the formula for a geometric sum. 

Moreover, this convergence is uniform on any closed sub-disk of |A| < 27r, because L(A) is a power 
series and e^^dx is bounded away from and from oo. Recall that differentiation commutes with 
uniform limits of holomorphic functions (as a consequence of the Cauchy formula). It follows that 
the derivative ^ commutes with the infinite order differential operators L(^-) on the right hand 
side of ©. Comparing the Taylor coefficients of A" on the left and right hand sides of @, we get 
a similar formula for f{x) = x", and hence for all polynomials. 

Remark 7. In higher dimensions, the problem of obtaining a formula for polynomials from a formula 
for exponentials is addressed in |KP2j . 

Weighted sums and ordinary sums. The Todd function is defined by 

Td(5) = — 

1 — e ^ 

The Bernoulli numbers are the coefficients in 

Td(-5) = ^ = 1 + i^.s^ 

k>l 

(We are following the conventions in Bourbaki Bo .) Since 

S/2 , -S/2 I / S S 

= i(Td{S)+Td(-S)), 

the coefficients b2k in the power series expansion (jSJ are the even Bernoulli numbers. Since 

Td(5) - Tdi-S) = — ^ - = S, 

1 — e ^ e"^ — 1 

we have 62^+1 = for all /c > 1, so the only difference between Td(5') and L(S') is the absence of the 
linear term in L(S'). Replacing L(-) by Td(-) on the right hand side of the Euler Maclaurin formula 
results in a formula for the ordinary sum of the values of / over the integers in [a, b]. However, 
we work with the formula for the weighted sum because this formula avoids "boundary effects" 
which occur in the formulas for the ordinary sum. 



4 Y. KARSHON, S. STERNBERG, AND J. WEITSMAN 

Euler Maclaurin with remainder. As before, we will denote the truncation of the power series 
L(S') at the even integer 2k by 1?^{S). Then 



d 15^ 1 d"^^ 



is a differential operator with constant coefficients involving only even order derivatives. In partic- 
ular, if g[h) is a function with 2k continuous derivatives, then 



(9) l?\^)(9{h)) 



h=0 



h=0 

The classical Euler Maclaurin summation formula with remainder can be formulated in the 
following way. 

Proposition 10 (Euler Maclaurin with remainder for intervals). Let f{x) be a function with m > 1 
continuous derivatives and let k = [m/2\. Then 

(11) ^7 = L2'=(A)L2fc(^) r^''/(x)dx + {-ir~^ t PMf(^\x)dx 

with 

(12) Pm{x) = 

ml 

where Bm{x) is the mth Bernoulli polynomial (see below) and where {x} = x — \_x\ is the fractional 
part of X. Moreover, the function Pmix) is given by 

°° 2 cos 2'Knx 



(13) P..(.) = (-1)'-'E ,2,„p 

if m = 2k is even and by 

A\ n / N / 2sin27rnx 

(14) P2k+i{x) = E (2^n)2fe+i 

if m = 2k + 1 is odd. 

Up to minor changes in notation, this result is formula 298 in jKnj . 

Note that if / is a polynomial then ()11() becomes an exact formula when m is greater than the 
degree of /. 

Let us recall the proof of Proposition 1101 Consider the difference between integration and sum- 
mation as the difference between two distributions: 



Po{x) :=l-Y^6{x -k). 



Integrating this, with the choice of constant of integration such that the integral from to 1 of the 
result vanishes, gives a distribution Pi which is given by the function 

(15) p,(:r) = x-lx\-^ 

at non- integral points. This is the famous zigzag function studied by Gibbs Idij . Notice that Pi is 
1-periodic and odd. 

If / is a continuously differentiable function on [0, 1] , we get 



[' Pi{x)f'ix)dx= Pi{x)f{x] 
Jo 



f{x)dx = IfiO) + J/(l) - /' f{x)dx, 

Jo ^ ^ JO 
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SO 

-1 



i/(0) + i/(l) = / 



f{x)dx+ / Pi{x)f'{x)dx. 



Summing the corresponding expressions over the intervals [a, a + 1], [a + 1, a + 2], . . ., [b — l,b], 
where a < b are integers, gives 

(16) J] 7= f f{x)dx+ f Pi{x)f{x)dx. 

[a,b] 

This is the first of a series of expressions relating the sum to an integral. 

We now take successive anti-derivatives: for m > 2, define Pmix) inductively by the conditions 
that -^Pmix) = Pm-i{x) and Jq Pm{x)dx = 0. Then Pm{x) is 1-periodic and satisfies Pm{—x) = 
{—l)'^Pm{x). Also, Pm{x) is Continuously differentiable m — 2 times; in particular, it is continuous 
if m > 2. Starting from (|16)) . we integrate by parts: 



[a,b] 



Y.'f = f f{x)dx+ P2f'\l- f P2{x)f"{x)dx 

= f f{x)dx+P2f\l-Pd"t+ f Pz{x)f^^\x)d 
J a J a 



X 

b 



^ f{x)dx+ P2ft - P3f"\l+ Pd^^'^l- 1^ PA{x)f^'Kx)d 



Noting that Pn{a) = Pn{b) = Pn{0) and i-2n+i(0) = (because P2n+i is odd), we get, setting 

k = [m/2\, 



(17) 5^7=/ fix)dx + P2{0)ft + Pmf^'^t + --- + P2k{0)f^'''''>t 

la,b] 

+ i-ir-' f Pra{x)f^'^\x)dx. 
J a 

Consider the polynomial 

(18) L[2'=1(5) := 1 + P2(0)52 + P4(0)S^ + . . . + P2fc(0)52'=. 
From (|17j) we get the remainder formula 

a a rb±h2 rb 

(19) 5^7 = L[2^](^)L[2^](^) / f{x)dx + (-1)—^ / Pm{x)f^'^\x)dx 

1^ dhi dh2 Ja±h, rn=h2=o J'^ 

for a function / of type C™", where k = \rn/2\. 

This formula becomes exact when / is a polynomial and when m is sufficiently large. We therefore 
have, by comparison with Equation @, 

(20) lI^'^I = l?^. 

This and ^ give (fTT|) . 

It remains to derive the expressions (fT2]). (|T3|) . and ([Ti|) for the functions Pm{x). 
The Bernoulli polynomials Bm{x) are characterized by the properties 

d 

Bo(x) = 1 , —Bmix) = mBm-i(x) , and / Bm{x)dx = for m > 1. 
dx Jo 



6 Y. KARSHON, S. STERNBERG, AND J. WEITSMAN 

In particular, 

Bi{x) = X - ^ = Pi{x) for all < x < 1. 

(See (|T5|) .) Integrating, we get that Pmix) = Brn{x)/m\ for all < x < 1, and since Pm is 
1-periodic, we get (fT^ . 

The Poisson summation formula says that, as distributions, 

fcez nez 

so 

oo 

Po{x) = 1 — 5{x — k) = —2 cos 2'Knx, 

keZ n=l 

as distributions. Integrating this, with the choice of constant of integration such that the integral 
from to 1 of the result vanishes, gives the expressions ()13() and (|14() . 

The idea of using these Fourier expansions goes back to Wirtinger ,W]. See Knopp |Knj . pp. 521- 
524. 

Remark 21. From H18|) and (|13() we see that the coefficients of the polynomials ^'^'^1(5') are 

where is the Riemann zeta function. Comparing this with the coefficients of the polynomials 
L'^'^iS), we get, from (j2ni) and ©, that 

reproducing Euler's famous evaluation of the Bernoulli numbers in terms of the Riemann zeta 
function. 

Proposition ll()| when applied to a C™ function f{x) of compact support, implies a similar formula 
for an infinite ray: 

(22) i/(a) + /(a + l) + /(a + 2) + ... 



a /•oo />oo 

/ fi^)d^ +(-ir-M P™(x)/(")(x)dx 

On Ja±h h=0 Ja 



for either choice of it, where k = [m/2j. Indeed, we need only choose b so large that the support 
of f{x) is contained in the set {x < b}, and then apply Proposition llOL Conversely, if we know 
((221) ^'^^ functions of compact support, then we can conclude (|11|) . Indeed, multiply / by a smooth 
function of compact support which is equal to one in a neighborhood of [a, b] and observe that 

p3) E'=E'-E' 

[a,b] [a,oo) [b,oo) 

when applied to a function of compact support. 

Twisted Euler Maclaurin with remainder for a ray. In extending the Euler Maclaurin for- 
mula to higher dimensions, we will need an expression for the "twisted weighted sum" 



oo 

-/(0) + j;A"/(r 



n=l 

when A ^ 1 is a root of unity, say, of order N, in terms of the integrals of /. 
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Let 



(24) 



Qo,a(x) = -^A"J(x-n). 



This is an iV-periodic distribution since = 1. We will take successive anti-derivatives of this 
distribution, with the constants chosen so that the integrals from to vanish: 



(25) 



-^Qm,x{x) = Qm^i,x{x) and J Qm,x{x)dx = 0. 



With this choice of constants, Qm,x{x) is A^-periodic for each m. 
Let us now look more closely at Qi^x{x). We have 



d 

dx 



l[n,n+i] {x) = 5{x - n) - 6{x - {n + 1)), 



so 



x-n) = ^—^Qo,x{x). 



Thus, 



n+l]( 



because the right hand side is periodic of period A'^ and its integral over [0, A'^] vanishes since 
1 + A + + • • • + A'^~^ = 0. If / is a continuously differentiable function of compact support, we 
have 



Qi,x{^)f'{x)dx 



A 



n=0 



n.+l 



A 



1-A 



/(O) + A/(l) + A2/(2) + 



so 



(26) 



l/(0)+J]AV(n) 



n>l 



I + -A^ 1 /(O) + Qi,x{x)nx)dx. 



2 1-A 



Successively applying integration by parts to H26() . we get 



1 



/(0) + A/(l) + AV(2) + 



^ + 1^1 /(O) - ^?2,a(0)/'(0) 



Q2,x{x)f"{x)dx 







^ + -r^ ] /(O) - g2,A(0)/'(0) + Q3,a(0)/"(0) + Q3,x{x)f^'\x)dx 



2 1-A 



I + /(O) - Q2,a(0)/'(0) + . . . + (-l)'=-iQ,,,(0)/('^-i)(0) 



+ -1 



Since 



(_l)"^-i/{"^-i)(0) 



5/i 



m /-co 



f{x)dx 



h=0 



we have thus proved this "twisted Euler Maclaurin formula for a ray" : 
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Proposition 27. Let 

M'''\S) = Q + S + Q2,x{0)S^ + Q3,x{0)S' + ■■■ + Qk,x{0)S\ 

for a root of unity A 7^ 1 . Then 

1 8 r°° 

(28) -/(O) + A/(l) + AV(2) + ■ • • = M'^'^^) J ^ f{x)dx 



h=0 

/•oo 

+ {-1)"-' / Qk,x{x)f^'\x)dx 
Jo 

iffeC^iR). 

Remark 29. Another twisted Euler Maclaurin formula for a ray appeared in jTj. 

We will now show that the polynomials that appear in Proposition EZI have the following sym- 
metry property: 

(30) M"^'^"'(5) =M'"'^(-5). 

The linear term transforms according to (|3U() because 



2 l-A-i 2 V 1-Ay \2 1-A 

For the other terms to transform correctly, we need to check that 

(31) Qm,x-^iO) = {-irQm,x{0) for all m > 2. 

As in the non- twisted formula, it is convenient to work with the Fourier expansions of the Qm,A's- 
Let 

A = e2-^. 

We are assuming that A ^ 1, and so j ^ mod A^. Then we can write the distribution Qo^x{x) as 



[X — n) 



Writing 

we see that the Fourier series of Qq^x{x) is 

Qo,A(x) = -5]e2-(^+'-)^ 

The indefinite integral, chosen so that the integral over [0, A^] vanishes, is obtained by dividing each 
Fourier summand by the coefficient of the exponent, so we have the following Fourier series: 



(32) Q^^^{x) = -Y, 



rez y2m{jj+r] 

Setting X = 0, replacing j by —j, and replacing r by —r in the sum, pulls out a factor of (— 1) 
which gives the required equation ()31() . and which implies H30|) . 
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Remark 33. For A = 1, if we define 

MM(5) = L2Lfe/2j(^) and Q,,i = Pk, 

then (|28() boils down to (|22() . So, with this notation, ()28j) also holds for A = 1. Notice that if A 7^ 1 
then M'='^(5') is a multiple of S, and that if A = 1 then M^'^(5) = 1+ a multiple of S. Finally, the 
symmetry property H3U|) continues to hold for A = 1 because the polynomials h"^^ are symmetric. 

3. The polar decomposition of a simple polytope 

In this section we decompose a polytope into an "alternating sum of polyhedral cones" . See 
or We will give a "weighted version" of this decomposition. 

Our polytopes are always compact and convex. A compact convex polytope A in M" is a compact 
set which can be obtained as the intersection of finitely many half-spaces, say, 

(34) A = Hin...r\Hd. 

We assume that (jSU is an intersection with the smallest possible d, so that the Hi's are uniquely 
determined up to permutation. We order them arbitrarily. The facets (codimension one faces) of 
A are 

ai = AndHi , i = l,...,d. 

Alternatively, a compact convex polytope is the convex hull of a finite set of points in M"'. If we 
take this set to be minimal, it is uniquely determined, and its elements are the vertices of A. 
For each vertex u of A, let 1^ C {1, . . . , d} encode the set of facets that contain v, so that 

i ^ ly if and only if v £ ai. 

Assume that A is simple, so that each vertex is the intersection of exactly n facets. For each 
i € ly, there exists a unique edge at v which does not belong to the facet af, choose any vector 
ai^y in the direction of this edge. (At the moment, these "edge vectors" are only determined up to 
a positive scalar. Later, when discussing integral polytopes, we will make a specific choice of the 
edge vectors.) 

The tangent cone to A at u is 

Cy = {v + r{x - v) I r >0 , X e A} = V M>oaj,^. 

We will "polarize" these tangent cones by flipping some of their edges so that they all "point in the 
same direction". This direction is specified by the choice of a "polarizing vector": a vector ^ G R"*, 
such that {^,aj^y) is non-zero for all v and j. With this choice, we define the polarized edge vectors 
to be 

tt _ J 0(i^y if {^,ai^y) < 0, 



(35) 

and the polarized tangent cone to be 



-tti^y if (C,aj,i,) > 0, 



(36) Cl = v+ 

We define the "weighted characteristic function" , 

(37) lX(x), 

to be the function on M" that takes the value on the exterior of A, the value 1 on the interior 
of A, and the value 1/2^^ on the relative interior of a codimension k face of A. So, for example, 
for an interval [a, 6] on the line, the function V^^^{x) assigns the value one to points a < x < b, 

zero to points outside the interval, and ^ to the points a and h. We use a similar definition for a 
polyhedral cone. 
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A =1/4 

• =1/2 

■ =1 

Figure 1. Polar decomposition of a triangle 

Proposition 38 (Weighted polar decomposition of a simple polytope). Let A be a simple polytope. 
For any choice of polarizing vector, we have 

(39) lX(x) = ^(-l)#n-,(x), 

V 

where the sum is over the vertices v of A, where Cf, is the polarized tangent cone, where IX(^) 0'''^d 
I'^j (x) are the weighted characteristic functions, and where denotes the number of edge vectors 

at V whose signs were changed by the polarizing process (|!-{5|) . 
This theorem is illustrated for the case of a triangle in Figure ^ 

Proof. We will prove this equality in two steps. First, for each x we will find a polarization such 
that the equality holds. Second, we will show that the right hand side of the equality is 
independent of the choice of polarization. 

Suppose that x A. Let ^ be a polarizing vector such that x) > y) for all y £ A. Then 
all the polarized cones "point away from x" and do not contain x. Formula (|39() for the polarizing 
vector ^, when evaluated at x, states that = 0. 

For every polarization there exists exactly one vertex v for which Cj, = C^,, namely, the vertex 
V such that v) is maximal. Conversely, for every vertex v there exists a polarization such that 

Cy = Cy. 

Suppose that x € A is in the relative interior of a face F of codimension k. Let v be any vertex 
of F. Let ^ be a polarizing vector such that Cf, = C^,. Let Fy be the codimension k face of Cy 
that contains F. Then l^(x) = I'^j (x) = . For each other vertex, v', the cone C^, is disjoint 
from the relative interior of F, so l"'^ (x) = 0. Formula (|39() for this polarization, when evaluated 

at X, states that (|) = (^) . 

To keep track of the different possible choices of polarization, we let 

El, ... , En 
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denote all the different codimension one subspaces of M"* that are equal to 

af^^ = {r]£ M"* I {r,, a^,^) = 0} 

for some j and v. (For instance, if no two edges of A are parallel, then the number of such 
hyperplanes is equal to the number of edges of A.) A vector ^ can be taken to be a "polarizing 
vector" if and only if it does not belong to any Ej. The "polarized cones" Cf, only depend on the 
connected component of the complement 

(40) W*\{EiU...UEn) 

in which ^ lies. Any two polarizing vectors can be connected by a path in R"* which crosses the 
"walls" Ej one at a time. We finish by showing that the right hand side of formula does not 
change when the polarizing vector crosses a single wall, Ek- 

As crosses the wall Ek, the sign of the pairing {^t^oij^v) flips exactly if E^ = aj-^^. For each 
vertex v, denote by S^(x) and S^(x) its contributions to the right hand side of formula before 
and after crossed the wall. The vertices for which these contributions differ are exactly those 
vertices that lie on edges e of A which are perpendicular to Ek- They come in pairs because each 
edge has two endpoints. 

Let us concentrate on one such an edge, e, with endpoints, say, u and v. Let denote an 
edge vector at v that points from f to u along e. Suppose that the pairing {^t,cte) flips its sign 
from negative to positive as crosses the wall; (otherwise we switch the roles of v and u). The 
"polarized tangent cones" to A at v, before and after crosses the wall, are 

Ci = v+Y^ R>oal^ + M>oae and 
(C« )' = v+J2 ^>o(^iv - I^>oae 

where le C {1, . . . ,d} encodes the facets that contain e. (The a^-^ are the same for the different 
^('s because the pairings {(,t,ct^jv) do not flip sign when crosses the wall for j G If,-) The cones 
cf, and (Cf,)' have a common facet and their union is 

Cl~v + Y, K>oaS> + Mae. 

This union only depends on the edge e and not on the endpoint v. (This uses the assumption that 
the polytope A is simple and follows from the fact that aj^u £ ^>o(Xj,v + Mog.) 

The contributions of v to the right hand side of (fS^ before and after crosses the wall are 

(41) S4x)=£l^„ and S;(x) = -elj;^„^, 

where e € {—1, 1}. Their difference is plus/minus the weighted characteristic function of C|: 

(42) S,(x)-SUx)=el^, 

The contributions of the other endpoint, n, have opposite signs than the respective contributions 
(|4H) of V, and their difference is minus/plus the characteristic function of C|. Hence, the differences 
Sv{x) — S^(x) and St((x) — S(^(x), for the two endpoints u and v of e, sum to zero. □ 

4. EULER MACLAURIN formula with REMAINDER FOR REGULAR POLYTOPES. 

We extend our notation (0) to a weighted sum over a simple integral polytope A: 
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where is the weighted characteristic function, introduced in 1)3 7() . which is equal to 1/2^ 

when X lies in the relative interior of a face of A of codimension k. We use similar notation for the 
weighted sum of a compactly supported function / over a simple polyhedral cone. 

With this notation, Proposition 1381 gives the following decomposition. Let ^ G M"* be a "polar- 
izing vector" . Then 

(43) ^ 7 = E(-i)*' E '/ 

where we sum over the vertices v of A, where Cf, is the polarized tangent cone to A at (see H36|) l. 
and where ij^v is the number of edge vectors at v that are flipped by the polarization process (|35() . 

For the standard closed orthant O = nr=i^>o ™ ^^"i have 

E'5= E E '^^("^1' 



rrir, 



onz" 



for any function g of compact support. 

By performing n iterations of (|22() . we obtain an Euler Maclaurin formula with remainder for 
the standard orthant: Let m be an integer and let k = [m/2j. If is a (7™" function of compact 
support and m > 1, then for any choice of the ±'s we have 



(44) 



n 



onz" 



i=l 



d_. 



g{x)dx 



0(±/ll,...,±h„) hl = - = hr,=0 

where 

O{hi,...,hn) = {t\ti + hi>0yi} 
denotes the shifted orthant, and where the remainder term is given by 

(45) R'^ig)^ 



/C{l,...,n} 



2k t 



0(±/ii,..,±?t„) 



40/ 



_d_ 

dxi 



g{x)dxi ■ ■ ■ dx,r, 



h=0 



(46) 



-dxi ■ ■ ■ dXn- 



This remainder can also be expressed as a sum of integrals over the orthant of bounded periodic 
functions times various partial derivatives of / of order no less than m and no more than mn. This 
fact follows from the formula 

— / (fdXi ■■■dXn = - / ^ 

ohi J J oxi 

0{hl,...,h„) 0{hl,...,hn) 

We can apply iterations of this formula with i G / to the /'th summand of H45|) because the 
non-smooth functions Pm are only applied to the variables xj for j ^ I. We get 

(47) R^M= E (-i)^"^-'^^"-i'i^/nL''(-^)n^-(^^)n(^™^)^^i---dx„. 



/C{l,...,n} 



O 



i&I 



10 



A regular integral orthant C is the image of the standard orthant O via an affine transformation 
of the form 



(48) 



,tn) ^ X = V + tiai + . . . + tnttn 
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where ai, . . . , a„ generate Z" and where v € Z". If ui, . . . ,Un & M"* is the dual basis to ai, . . . ,an 
then the image of 0(/ii, . . . ,hn) under this transformation is given by the inequahties 

(49) {ui,x)-{ui,v) + hi>0. 

We denote this expanded orthant by C{h). If / is a C""" function of compact support and 

g{ti, ...,tn) = fihai + . . . + tnOn) 
is its puhback under the transformation H48|). then 

^ / = ^ 5 and / f{x)dx = / g{t)dt, 
and so we have an Euler Maclaurin formula for regular orthants: 



(50) J2'f=Il^''(^) f f(^)dx 



d_ 

dhj 

C{±hi,...,±h„) 



hi=---=hn=0 



where 



Let A be a compact convex polytope. As in Section |31 A can be written as an intersection of 
half-spaces 

(51) A = n . . . n Hd, where Hi = {x \ {ui,x) + > 0} for i = 1, . . . , d 

and d is the number of facets of A. The vector Ui € M"* can be thought of as the inward normal 
to the ith facet of A; a-priori it is determined up to multiplication by a positive number. If all the 
vertices of A are integral, then the Uj's can be chosen to belong to the dual lattice Z"*, and we 
can fix our choice of the Uj's by imposing the normalization condition that the Uj's be primitive 
lattice elements, that is, that no Ui can be expressed as a multiple of a lattice element by an integer 
greater than one. (The fact that a normal vector u to a facet a can be chosen to be integral is a 
consequence of Cramer's rule. Indeed, we can choose integral edge vectors . . . , /3„ that emanate 
from a vertex on a such that Pi, ■ ■ ■ Pn-i span the tangent plane to a and is transverse to a. 
Solving the linear equations {u, (3i) = ... = («, (3n-i) = and {u, fin) = 1, we get an inward normal 
vector u with rational entries; clearing denominators, we may assume that u is actually integral.) 

We can then consider the "dilated polytope" A(/ii, . . . , hd), which is obtained by shifting the ith 
facet outward by a "distance" hi. More precisely, 

d 

A(/i) = I (uj, x) + /ii + /ij > 0} where h = {hi, . . . , h^). 

4 = 1 



Now assume that A is simple. Then A{h) is simple if h is sufficiently small. The polar decom- 
position of A(/i) involves "dilated orthants". However, dilating the facets of A outward results 
in dilating some facets of C| inward and some outward. Explicitly, for i ^ ly — {h? • • • ^^n 

}, the 

inward normal vector to the ith facet of Cf, is 



(52) u. 



l,V 




if al^^ = Oi^y 
if 4,v = 
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Hence, the dilated orthants that occur on the right hand side of the polar decomposition of A(/i) 
are c!,(/if^ . . . , /if^^J, where 



(53) 



hi 



if O^j ^. — Oii y 



-hi if a - , 



-OLi, 



This subtlety in the signs does not effect the formula for regular polytopes, because of the symmetry 
of \?^\ however, it is needed to derive the formula for simple polytopes. 

If A is a regular integral polytope, then the C|'s are regular integral orthants, and the "dilated 
orthants" are exactly as in (|49|) . We then have, by (|4.S|) . (|50j) . and the symmetry of L^*^, 



AnZ" 



E(-i; 



n 

«e/t,={ii,...,i,i} 



f{x)dx 



cU±h,^,...,±h,„) 



\ 



h=0 



J 



We may multiply the differential operator O-'-'^'^C^) ™ above expression by any number of 
operators of the form L^'^(g^) where j ly, since all that will remain of these operators is the 
constant term 1, all actual differentiations yielding zero. The right hand side of H54() is then equal 
to 



nL-,A)E(-i)*' / 



f{x)dx 



Cl(h*. ) 



h=0 



nL-(^) / 

1 = 1 A /I. I. \ 



A{hu...,ha) 



h=0 



where 
(55) 



S'Sif) :=E(-1)*'^™ (/) 



and where {ii, . . . , in} = Iv in the u 'th summand. 
Notice that both 



E'/ 



AnZ" 



and 



A(/ii,...,/id) 



i=l 



h=0 



do not depend on the choice of polarization, and both vanish on any function / whose support 
is disjoint from the polytope. So the same must be true of the remainder. We have proved the 
following result: 
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Theorem 1 f |KSWj l. Let m > 1 be an integer. Let A C M" be a regular integral polytope and f 
a C""" function of compact support on M". Choose a "polarizing vector" for A. Then 



AnZ" 



f{x)dx 



A{hi,...,ha) 



+ S^{f) 



h=0 



where k = [m/2j and where S'^{f) is given by (|55|1 . This remainder can be expressed as a sum 
of integrals over orthants of bounded periodic functions times various partial derivatives of f of 
order no less than m and no more than mn. Finally, this remainder is independent of the choice 
of polarization and is a distribution supported on the polytope A. 



5. Finite groups associated to a simple integral polytope and its faces. 

In order to extend Theorem ^ to simple integral polytopes that may not be regular, we must 
extend the Euler Maclaurin formula for regular orthants (|50|) to a formula that is valid for simple 
orthants that may not be regular. In this section we analyze certain finite groups that arise in this 
generalization. 

Let C be a simple integral orthant. This means that we can write C as the intersection of n 
half-planes in general position, 

(56) C = /fi n . . . n Hn where Hi = {x \ {ui, x) + fii > 0} for z = 1, . . . , n, 

and that the vertex of C is in Z". The Uj's are inward normals to the facets of C, and we choose 
them to be primitive elements of the dual lattice Z"*. (See Section We choose ai, . . . , a-n to be 
the dual basis, that is, 

/ \ j'^ j = ^ 

{uj,ai} = < 

The aj's are edge vectors of C, but they might not be integral: they generate a lattice in M" which 
is a finite extension of Z". This extension is trivial exactly if A is regular at v, that is, if the Uj's, 
generate the dual lattice Z"*. To the cone C we associate the finite group 

(57) r:=Z"7^Zni. 
So this group is trivial exactly if the cone C is regular. 
Lemma 58. Ln the above setting, 

(59) 7 ^ e2^'<^'^> 

is a well defined character on T whenever x = miui + . . . + ninCtn where mj are integers. This 
character is trivial if and only if x ^17^ . 

Proof. Let 7 G Z"* be an element that represents 7. If we expand it as a combination of the basis 
elements Uj, so that 

7 = biui + . . . + 6„n„, 

then (7, aj) = bj for all j. If 7' is another element of Z"* that represents 7, then it differs from 7 
by an integral combination of the Uj's. So {'y',aj) differs from (7,aj) by an integer. Hence, when 
x is an integral combination of the CKj's, the pairing (7, x) is well defined modulo Z, and so e'^'^''^"''^'^ 
is well defined. 

Finally, the character H59|) is trivial if and only if (7, x) is an integer for all 7 € Z"*, and this 
holds if and only if 2; S Z". □ 



16 Y. KARSHON, S. STERNBERG, AND J. WEITSMAN 

Let A be a simple integral polytope in M", given by ()51|) . For each vertex v of A, let C 
{1, . . . ,d} denote the set of facets of A that meet at v. The normal vectors Ui, i £ I^, form a basis 
of M"*. We choose 

(60) Oi,^ G M" , ielv 

to be the dual basis. 

Given any face F of the polytope A, let Ip denote the set of facets of A which meet at F. 
Because A is simple, the vectors Ui, for i S Ip, are linearly independent. Let Np C M"* be the 
subspace 

Np = spanjuj | i € Ip}. 

Remark 61. It is natural to define the tangent space to the face F to be TF = span{x— y \ x,y G F}, 
the normal space to F to be the quotient W^/TF, and the co-normal space to be the dual of the 
normal space. With these definitions, Np is the co-normal space to F. 

To each face F of A we associate a finite abelian group Tp. Explicitly, the lattice 

Vp=Y^ Zui C Np 

ieip 

is a sublattice of Np n Z"* of finite index, and the finite abelian group associated to the face F is 
the quotient 

(62) Tp := {NpnZ''*)/Vp. 

If F = f is a vertex, this is the same as the finite abelian group associated to the tangent cone C„ 
as in (|57|) . 

Let E and F be two faces of A with F Q E. This inclusion implies that Ip ^ Ip, and hence 

(63) {ui}i(zi^ C {ui}i(zip. 

Because these sets are bases of the vector spaces Np and Np, we have an inclusion 

Np C Np. 

Because the sets occurring in are Z-bases of the lattices Vp and Vp, we have A''^; nVp = Vp. 
Hence, the natural map from Tp = (Z"* n Np)/Vp to Tp = (Z""* D Np)/Vp is one to one, and it 
provides us with a natural inclusion map: 

if F C F then Tp C F^. 

We define a subset F^ of F^? by 

(64) F^:=F,.\ J Tp. 

faces E such that EDF 

Then 

(65) F, = y F^. 

{F:veF} 

Recall that 

(66) A^j> := e^^^^T'^J."), for 7 G F^ and j G /„, 

is well defined, by Lemma EHl It is a root of unity. We will need the following results. 
Claim 67. //7 G Tp and j G Ip, then X^yj^v is the same for all v (z F. 
This allows us to define X^,j,p for 7 G F^? and j G Ip such that 

^7,j,P ~ for 7 G Fi? and j G Ip, if v £ F. 

Claim 68. If ^ & Tp and j G Iv~^If then X-yj^v is equal to one. 
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This allows us to define X^j^p = 1 when 7 G F^? and j G {1, . . . , d}\Ip. Then 

(69) '^'y,j,F = X'y,j,v for 7 G Fi? and I < j < d, if v € F 
and 

(70) Kj,F = l ioT^eTpiij ^If. 
Claim 71. If j E F^ and j ^ If, then X^yj^F 7^ 1- 

Proof of Claims [F^ and\71\ Let 7 G Tp be represented by 

7 = ^ biUi G tVf n z"* 

for some bi G M. (See 

Let V G F. Because aj^^, j G It,, is a dual basis to Uj, j £ 1^, we have 



(7, Oj.t;) = j 

Hence, 



6j j G /f 

j G Iv^If- 

~ I 1 i G 



is independent of v and is equal to 1 if j G ly^Ip. This prove Claims EZI and EHl 
Let j £ Ip. If A-yj^i? := e^'^*^^ is equal to one, then bj is an integer, so 

(72) 7'= 6,ni 

also represents 7. Let -E D F be the face such that Ip = Ip^lJ}. Then, by (f7^. 7 G F^;. This 
proves Claim [7T1 □ 



6. EULER MACLAURIN formula with REMAINDER FOR SIMPLE INTEGRAL POLYTOPES. 

Let us begin by deriving an Euler Maclaurin formula with remainder for a simple integral orthant. 

We recall the set-up of Sectional Let C be a simple integral orthant. Let ui, . . . ,Unhe the inward 
normals to its facets, chosen to be primitive elements of the dual lattice Z"*, and let ai, . . . , be 
the dual basis to the Uj's, so that 

n 

i=i 

Let 

F = Z"7^Zuj 

be the finite group associated to C. By Lemma 1551 7 1— > e'^'^^^"''^'^ defines a character on F whenever 
X G ^Zoj, and this character is trivial if and only if x G Z". By a theorem of Probenius, the 
average value of a character on a finite group is zero if the character is non-trivial and one if the 
character is trivial. So 

|F[ ^ |o ifx0Z" 
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for all X € ^Zqj. Then, for a compactly supported function f{x) on M", 

cnz" X y I I -yer J 

(73) = ^^^'e2-<^'-)/(x) 

where we sum over all 

(74) X = V + miai + . . . + mnCtn, 

where the mj's are non-negative integers. 

The cone C is the image of the standard orthant O under the affine map 

(75) {ti, . . . ,tn) ^ X = V + tiai + . . . + tnOn- 

This map sends the lattice Z" onto the lattice Yl '^cnj ■ The inverse transformation is given by 

(76) ti = {ui,x- v) . 

Let us concentrate on one element 7 G F. Because u E Z", from ()74|) we get 

n 

i=i 

so that the inner sum in 1)73^ becomes 

(77) E V-- - E'^n"5("^i'---'"^«)' 

mi>0 m„>0 

where 

g{ti, . . . ,tn) = f{v + hai + . . . + tnan). 
Recall that we had the twisted remainder formula 

POO 

h=0 JO 



^'A™9(m) = M'^'^(|-) Hgm 
m>0 ■^^'^ 



for all compactly supported functions g{x) of type C^, where k > 1, where A is a root of unity, and 
where M^''^ is a polynomial of degree < k. (See p8|) and Remark 1331 ) 
Iterating this formula, the sum in (|77|) can be written as 

^"'^'(7.7-)/ •••M'='^"(— ) / giti,...,tn)dh---dtn + Rl\Xu...,Xn;g) 

ohi J^h-, dhn J-h„ 

(78) =n^*^'^'(a^) J 9iti,---,tn)dti---dtn + Ri\Xi,...,Xn;g) 

0{/i) 

with 

0{h) = {{ti, ...,tn) \U> -hi for all i} 
and where the remainder is given by 

i?fc*(Ai, . . . ,Xn;g) = 

Qk 



/c{i,...,n} ie/ * o;^^) i^/ iii » 



/i=0 
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Using (|46() to express this as a sum of integrals over the (non-shifted) orthant, we get 
(79) Ri\Xi,...,Xn;g) = 



Qk 

7C{l,...,n} Q iel i^I 10 « 



^ j Y\ M^'^>(-|-) n Qk,xM n ^5(ti, . . . , in)'iii • • • 



We now perform in ((7SJ and 1)79(1 the change of variable given by the transformation (|76(1 . The 
integrals over 0{h) and O get replaced by integrals over C(/i) and C times the Jacobian factor |r|, 
the function g{t) gets replaced by the function /(x), and the partial derivative gets replaced by 
the directional derivative Dq-. Summing the expressions (|78|1 over 7 € T and dividing by |r|, we 
get, by JZSl), 



cnz" 7Gri=i 



h=0 



where 

and where 
(81) R^U):-- 



7er/C{l,...,n} Q iG/ 40/ j^/ 



Let A be a simple polytope, given by ((^ . Choose a polarizing vector for A and let Cf, denote 
the polarized tangent cones. The inward normals to the facets of CE are given by the dual 
basis to these vectors is aj^, i G /„, and the roots of unity that appear in the Euler Maclaurin 
formula for Cf, are then 

Also recall that the polar decomposition of A(/i) involves the dilated orthants Cl{h\^ v'' ■ ■ ■ ■> v) 
where h\ ^ are as in . 

Let A: > 1 be an integer. For any compactly supported function / on M" of type C^^ , the polar 
decomposition of A(/i) and the formula (|8U() give 



(83) E = E(-i)*'^ E 



AnZ" 



where the remainder is given by 

(84) Rtif):=Y^{-l)*^R'^hf). 



f{x)dx 



h=0 
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Note that either hi = and A«,,,, = A,,,., or hi = -h, and X^^^ = A"],,. By the symmetry 
property (PU)) . this gives 

dhl' 'dhj 

For j ^ ly, because A^J,^, = 1 (see (IZOl)), we have M^^'^-''^-" {-^) = 1+powers of Also, the 
cone Cl{h\^ v' ■ ■ ■ ' v) independent of hj for j ^ ly Therefore, (|HH|) is further equal to 



(85) 



76r„ j=i 



where in the T;'th summand {zi, . . . , in} = ly 
Define 



(86) 



Then we have, by H69|l . 
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/i=0 



M^^p = Yl M^'^^.^.^( — ) for 7 G TiT 



(87) p = „ whenever 7 G and w G F, 

where we identify 7 G F^? with its image under the inclusion map F^? ^ F^ 
Then (jHSJ is equal to 



J](-l)#^ ^ M^,, / /(x)cix 



?i=0 



f{x)dx 



h=0 



by ()65() and 1)8 7|) . In the interior summation we may now add similar summands that correspond to 
V ^ F. These summands make a zero contribution to ()88p for the following reason. If v ^ F then 
there exists i G Ip^Iv Because i ^ 1^ = {ii, . . . , in}, the cone cS(/i^^ vi ■ ■ ■ independent of 

hi. So it is enough to show that 'M.^ p is a multiple of But because 7 G F^ and i £ Ip, we 
have X-y^i^F 1- (See Claim ITU) By Remark 1551 this implies that M'^''*'T'^'^(^), which is one of 
the factors in M^^, is a multiple of Hence, (|S5|) is equal to 



F ^ert aiii, 



f{x)dx 



h=0 



76r^ 



k 



f{x)dx 



h=0 



We have therefore proved our main result: 
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Theorem 2. Let IS. he a simple integral polytope in M" and let f G C" (M") be a compactly 
supported function on M" for k > 1. Choose a polarizing vector for A. Then 



h=0 



where M!^ ^ are differential operators defined in (|86j) and where the remainder R^{f) is given by 
equation \84^ . Moreover, the differential operators ^ are of order < k in each of the variables 
hi, . . . ,hii. Also, the remainder is a sum of integrals over orthants of bounded periodic functions 
times various partial derivatives of f of order no less than k and no more than kn. Finally, this 
remainder is independent of the choice of polarization and is a distribution supported on the polytope 
A. 

7. Estimates on the remainder and Euler Maclaurin formulas for symbols and for 

polynomials 

We first recall a definition from the theory of partial differential equations. A smooth function 
/ € C°°(M") is called a symbol of order N if for every n-tuple of non-negative integers a := 
(ai, . . . , o„), there exists a constant Ca such that 

|5?^..9^/(x)|<C7,(l + |x|)^-H 

where \a\ = 'Yli^i- particular, a polynomial of degree is a symbol of order N. Note that if 
/ is a symbol of order N on M" then its derivatives of order a are in if A < \a\ — n. In this 
section we will show that the Euler Maclaurin formula of Theorem |21 can be extended to symbols, 
and gives rise in this way to an exact Euler Maclaurin formula for polynomials. 

To make further progress, we first require an estimate on the remainder term R^{f). We recall 
(see Theorem|21) that this remainder can be expressed as a sum of integrals over orthants of bounded 
periodic functions time various partial derivatives of / of order no less than k and no more than 
kn. Explicitly, from (jHlJ and (|81() we get 

(90) = E(-i)*^ E 



/ Hm''"" )T[Q,^^ ((4.,^-^))n(^a« ?f{x)dx. 



ie/ j0 30 



The functions Qk,x are bounded and periodic, and M^'^{—Da.) are differential operators of order 
k. In particular, each integrand in the formula for R^{f) is dominated by a constant times 

sup \di'---di"f\ 

where the supremumis taken over all n-tuples {ji, . . . with k < + . . ■+jn < nk. Consequently, 
R^{f) is well defined (by the same formula 1^90^ ) for any smooth function / whose derivatives of 
order between k and nk are integrable on M". In particular, it is defined when / is a symbol of 
order less than k — n. Moreover, we get an estimate for the remainder: 

Proposition 91. The distribution R^{-) extends to symbols of order less than k — n, by the explicit 
expression ()9U|) above, and satisfies an estimate of the form 

\Rt{f)\ < K{k, A) • sup|,.^,...,,,^}|a^'^ . . . dtf\L,(R'^), 

where the supremum is taken over all n-tuples {ji, ■ ■ ■ , j„} with k < ji + ■ ■ ■ + jn 1^ nk. 

By applying this estimate we will obtain the following Euler Maclaurin formula for symbols. 
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Theorem 3. Let A be a simple integral polytope in R", let f be a symbol of order N on M", and 
choose k > + n + 1. Then 



(92) E '/ = E E ^',F [ /(^)^^ 

AnZ" F ^gr^ -'A(/i) 



h=0 



where are differential operators defined in (jHSJ and where the remainder term R^{f) is 

defined in (|Un|) . Moreover, the differential operators p are of order < k in each of the variables 
hi, . . . ,hd- Also, the remainder R^{f) satisfies the estimate of Proposition\91\ 

In the case where / is a polynomial, this formula gives rise to an exact Euler Maclaurin formula. 

Corollary 93. Let p be a polynomial on M", and choose k > deg p + n + 1. Then 



p = E E ^T.i^ / ^(^)^^ 

AnZ" F j(zri>^ JA(h) 



h=0 



Proof of Theorem\^- Let x be a smooth function on M" which is equal to one on some open ball 
about the origin that contains A and which is supported in some larger ball, say, of radius R. 
Define Xxi^) •= xi^/^) for all A > 1. Then the function 

fx ■= fXx 

is a smooth compactly supported function on M". We apply Theorem |2l to obtain 



(94) E /a = E E / hi^)d^ 



+ Rtifx 

h=0 



which is valid for any k > 1. 

Since fx equals / on a neighborhood of A if A > 1, the left hand side and the first of the two 
summands on the right hand side in ()94|) are equal to the corresponding terms in H92[). Thus, to 
deduce (|^^ from it suffices to prove the following claim: 

(95) Ihn R^ifx) = Rtif) if k>N + n + l. 

To prove this claim, we apply the estimate in Proposition |^ to the difference R^{f) — R^{fx) = 
R^{f{l — Xa))- We expand each of the derivatives appearing in the estimate into a finite sum of 
products of derivatives of / and derivatives of 1 — xa- The leading term has the form g{l — xx) 
where g is a derivative of / of order q with q > k. Because / is a symbol of order A^, the function 
g is dominated by a constant times the function x i— > (1 + Ixl)^""^. Because q > k > N + n + 1, the 
function g is in Li. It follows that the Li norm of g{l — xx) converges to zero as A ^ oo. Each of 
the remaining terms in the expansion has the form 

(96) A-^5(-)x(-/A) 

where 5 is a derivative of / of order q and x is a derivative of 1 — x of order s, and where s > 1 
and s + q > k. The function g is dominated by a constant times the function x 1— > (1 + |x|)^~'^; the 
function x is bounded; the function x('/A) is supported on the ball of radius XR about the origin. 
Hence, the Li norm of the product g{-)x{'/^) is bounded by a constant times X^-i+"- if — g + n 
is non-negative and by a constant otherwise. Since s > 1 and s-|-(7>A;>A^-|-n-|-l, the Li norm 
of the term (|96|) is bounded by some constant multiple of A~^. Letting A — > 00, we see that each of 
the terms in the estimate approaches zero as A — > 00. This implies (|95)) . Theorem |31 follows. □ 
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